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Abstract 

We compare the interaction parameters measured on LaMnOa to single site dynamical mean 
field estimates of the critical correlation strength needed to drive a Mott transition, finding that 
the total correlation strength (electron-electron plus electron-lattice) is very close to but slightly 
larger than the critical value, while if the electron lattice interaction is neglected the model is 
metallic. Our results emphasize the importance of additional physics including the buckling of the 
Mn-O-Mn bonds. 

PACS numbers: 71.10-w,71.30.+h,75.10.-b 
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The "colossal" magnetoresistance (CMR) materials are widely regarded as paradigm 
"strongly correlated" systems in which strong local interactions combine with orbital and 
lattice effects to produce a range of exotic behavior Many papers over the last decade 
have argued that LaMnOq, the "parent compound" of the CMR family, is a strongly cor- 
related Mott insulator 

mm 

, while others have argued that short or long ranged 
Jahn- Teller order is the key physics . The debate has continued in the literature up to 
the present. 

Two classes of difficulty have complicated the discussion of LaMnOs. The first concerns 
the definition of "Mott insulator" . While there is a general agreement that the term "Mott 



insulator" 



y, |7j refers to materials in which the electronic correlations are strong enough to 
lead to insulating behaviors in the absence of long ranged order or significant short ranged 
order, it is not easy to implement this definition in practice since most candidate Mott 
insulators exhibit some form of long ranged order at low temperature (in LaMnOs "low 
temperature" means below the orbital ordering temperature Too ~ 7507^) and neither long 
ranged order nor any short ranged correlation can easily be "turned off" experimentally. The 
second class of difficulty concerns the ambiguity of the concept of "correlation strength". 
The paradigmatic model for Mott insulator is the one orbital Hubbard model in which there 
is just one interaction parameter. The rich multiplet structure associated with partially filled 
d-shells in materials such as LaMnOs includes intra and inter orbital Coulomb repulsions, 
Hund's couplings and Jahn- Teller splittings. Understanding which interaction is the most 
important has not been easy. These difficulties can be circumvented theoretically. In a theory 
the effects of different interactions can be disentangled. Further the single site dynamical 
mean field theory (S-DMFT) |^ neglects intersite correlations entirely, but produces a 
metal-insulator transition for interactions larger than a critical value. We propose that a 
material is "strongly correlated" if, for the relevant interaction parameters, the single site 
DMFT approximation produces an insulating solution at zero temperature without long- 
ranged order. Materials with weaker but still non-negligible correlations should be referred 
to as having intermediate correlations. 

In this paper we combine experimental and theoretical information to show that the local 
interactions in LaMnOs are very close to the critical values needed to drive a metal-insulator 
transition in the single site DMFT approximation and in the absence of long ranged or short- 
ranged order. While aspects of our analysis have appeared in the literature, in this paper we 
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present a comprehensive view which allows us to resolve a dispute which has continued in 
the literature up to the present. The close proximity of the interactions to the critical value 
means that the behavior is extremely sensitive to the conduction bandwidth, to magnetic 
order, and to the electron-lattice coupling. In particular, for the bandwidth obtained from 
local density approximations for the observed structure the material is insulating in the 
paramagnetic phase but would be metallic in a hypothetically ferromagnetic phase or if 
the electron-phonon coupling were set to zero. Further, the observed crystal structure (at 
high temperatures above the orbital ordering temperature) differs from the ideal perovskite 
structure by a GdFeOs rotation. If this rotation is removed, the material becomes metallic. 

The rest of the paper is organized as follows. We first establish the Hamiltonian and 
present qualitative considerations. We next present an analysis of the changes in optical 
conductivity across the Neel transition which, following Ref 2] and RefQ], allows us to es- 
timate the interaction strength. We then present single site DMFT calculations for various 
model parameters which allow us to estimate the proximity of the parameters to the thresh- 
old values. We conclude with a summary and a general classification scheme. 

We now estimate the electronic parameters, beginning with the bandwidth W. LDA and 
LDA+U calculations reveal that the Cg bands are well described by a nearest neighbor tight 
binding model with hopping amplitude t jol; the Cg symmetry implies a directional structure 
to the hopping so that there are two bands and W = 6t. The bandwidth is found to depend 
strongly on the crystal structure; in particular on the magnitude of GdFeOa rotation away 
from the ideal perovskite structure The hopping increases from 0.5eV to 0.65eV when 
the GdFeOs rotation is decreased from the value observed for LaMnOa to zero (with the 
lattice parameters held constant). 

To estimate the interaction strengths we follow the analysis of Kovaleva et al \l\ who 
show that the optical spectrum of LaMnOs exhibits peaks which have a dependence on 
polarization and temperature allowing them to be associated to atomic-like excitations 61] 
with a reasonable degree of confidence. To establish notation and estimate uncertainties we 
give the details of the analysis here. The formulation and results are similar to those given 

n n 

by Kovaleva et al [2|], and the physical arguments were introduced in Ref [6|. 

We now turn to the local interactions. We assume (as is apparently the case in the 
actual materials) that the crystal field (ligand field) is large enough that the t2g levels are 
well separated from the Cg levels so that the pair hopping between t2g and Cg orbitals is 
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quenched (this assumption was not made in [2] which accounts for the differences between 
their results and ours). We also assume that the local interactions are strong enough that 
level electrons are in their maximum spin state, and may be treated as an electrically 
inert core spin of magnitude \Sc\ = 3/2. The on-site Hamiltonian in the manifold is then 

Hioc = '^{U - J)ni^^n2^^' + U'^ ^1,^^1,1 + J{c\^^c\^^C2,iC2,i + h.c.) 

a,a' i=l,2 

- 2 Jsi • 5*2 - 2JhSc ■ (si + S2) + A(ni - 722) (1) 

Here Si = J2af3 <^la^ai3Cii3, \Sc\ = 3/2 and A is a crystal field splitting of levels arising from 
a Jahn- Teller distortion of the Mn-Og octahedron which may be static or dynamic, and 
have long ranged order or not. For a free ion Jh = J', we assume this henceforth because 
the general consensus is that expect for U, intra atomic interactions are insensitive to solid 
state effect (screening). The eigenstates of Hioc are characterized by the particle number 
(n), total spin (Stot) and total Cg spin {Se^), and the orbital configuration (O). We label the 
two-electron states as '^''^!>+^0{Stot)- There are 16 1-electron and 24 2-electron eigenstates, 
taking the configurations of the core spin into account. 

In the atomic picture, in the ground state of LaMnOs each Mn atom is in the state 
n = 1, Slot = 2. The electron is in the particular orbital state picked out by the crystal 
field splitting. An optical transition then leaves one site in the state n = Stot = 3/2 and 
one site in the state n = 2 with orbital state labeled by O and spin state characterized by 
Stot = 5/2, 3/2, 1/2 and Se^ = 1, 0. Optical peaks are at energies AE{Stot, Seg, O) = E{n = 
2,Stot,S,g,0) + E{n = 0,Stot = 3/2) - 2E(n = 1,^^^ = 2, S^^ = 1/2). Table I lists the 
2-electron eigenstates, degeneracies and the corresponding optical transition energies. 



States 


Degeneracy 


AE 


^^2(5/2) 


6 


U - 3 J/2 + 2A 


^^2(3/2) 


4 


U + 7J/2 + 2A 


^^2(1/2) 


2 


U + 13 J/2 + 2A 


^^-(3/2) 


4 


U + 9 J/2 + 2A - V4A2 + J2 


1^(3/2) 


4 


U + 9 J/2 + 2A + V4A2 + J2 


^E+{3/2) 


4 


U + 7J/2 + 2A 



Table I: The 2-electron eigenstates of H (Eq([T])) la belled by Cg spin (superscript), 
representation (letter with or without subscript) and total spin degeneracy of the 
eigenstates (parenthesis), and the corresponding optical transition energies (see text) 
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The ^A2{5/2) are the states of maximal spin and are favored by Hund's rule. They 
necessarily have one electron in each orbital, hence a crystal field energy 2A higher than 
the starting state. The ^^2(8/2, 1/2) have the same orbital and Cg spin configuration as 
^^42(5/2) but lower total spin. The ^^42(1/2) state is not connected to the ground state by 
the optical matrix element and will not be considered further here. 

The ^£"+(3/2) state is the low eg spin configuration with one electron in each orbital and 
the remaining two states ^£'"(3/2) and ^^4(3/2) are Cg singlets made up of linear combina- 
tions of states with two electrons in the same orbital. These states are split by the crystal 
field but coupled by the inter-Cg pair- hopping. If 2 A » J, the pair hopping is quenched 
and we may identify the state ^E^ as the coming from the two electron state with both 
electrons in the orbital favored by the JT splitting and as the state with both electrons 
in the disfavored orbital. The state '^A is at high energy and is connected to the ground state 
by a very weak matrix element; it will be disregarded. The relevant portion of excitation 
spectrum of the exact model therefore consists of transitions to the high spin, '^E~ and ^E^ 
states implying peaks at 

AEhs = U-3J/2 + 2A 

AE^s = U + 9 J/2 + 2A- V4A2 + J2 

AE+g = U + 7 J/2 + 2 A (2) 

Note that the latter two peaks are degenerate at A = 0. 

An extensive theoretical literature exists on the problem of electrons coupled to classical 
core spins. We briefly discuss how to make comparison to these results. In the classical core- 
spin model we write the coupling between core spin Sc and conduction spin cxei {\(^ei\ = 1) as 
Hci = —JciSc ■ (^ei- The energy difference between high spin {ffei parallel to Sc) and low spin 
{aei anti-parallel to Sc) is 2JciSc. In the quantum model the energy differences depend on 
the total spin and electron number. The high-spin/low-spin difference measured in optics is 
5 J while the difference in energy between n = 1 high-spin/low-spin is 4 J. Therefore there 
is an approximately 20% uncertainty in the classical parameter 2JciSc- 

The critical U for a model with the interaction given in Eq([l]) has not been calculated. 
However there are several limits in which the manganite model maps on to an effective one 
orbital model; for these cases we may estimate the critical U . The first limit is of strong 
ferromagnetism. If the core spins are fully polarized the spin degree of freedom is quenched, 
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the orbital degree of freedom acts as a spin and the Jahn- Teller coupling A as a magnetic 
field. Projection of Eq([T]) on to the maximum spin manifold then yields a one orbital 
Hubbard model with a Uejf = U — 3 J/2 (if A = 0). For this model the critical U for the 
metal-insulator transition is Uc2 ~ 1.5W ~ 4.5eV [a [l^. A second simple limit is Jh ^ oo 
(the Jh in the actual materials is far from this limit). In this case, in the paramagnetic phase 
the spin degree of freedom is a gain quenched and Ueff is again U — 3J/2. The bandwidth 
is reduced by a factor of ^2 0,13] suggesting t/c2 ~ 1. 51^/^2 ~ I.IW ~ 3.3eV. 

Experiment on orbitally ordered LaMnOa |2| identifies two clear peaks at 2eV and 4.4eV, 
along with weaker features at 4.7eV and higher energies. In this material it is reasonable 
to regard the Jahn- Teller distortion A as frozen-in. For T < T/v ~ lAOK LaMnOa is an 
A-type antiferromagnet with ferromagnetic planes, which we take to define the x — y plane, 
antiferromagnetically alternating in the remaining, z, direction. Comparison of spectra 
taken with electric field along and perpendicular to z and at temperatures above and below 
T/v implies {2!, that the 2eV peak corresponds to the HS (high-spin) final state and the 
4.4eV to a LS (low-spin) final state. There is only one HS state so we identify U — 3 J/2 + 
2A ~ 2eV. There are two candidate LS states and therefore there is an uncertainty in the 
peak assignment. Consideration of the optical transition strengths implied by the observed 
orbital order suggests that the E^g state should be more prominent implying we identify 
2AeV=AE+g - AEhs = 5 J so that J = 0.48eV and f/ + 2A ~ 2.7eV. In this interpretation 
the transition to the AE~[g state would give rise to a weaker feature below the mean peak, not 
resolved as a separate excitation because of the band-broadening. The observed approximate 
0.5eV width then implies that the AEf^ cannot be more than about 0.5eV below AEfg 
implying 2A < O.TeV and 2.0eV< U <2.7eV. Alternatively we may identify the 4.4eV with 



Ej^g implying 2.4eV~ 6J — v4AM-J2. We would further identify the 4.7eV peak with Ej^g 
implying J ~ 0.54eV, 2 A ~ 0.64eV, and U ~ 2.18eV. The ambiguity in peak assignment 
therefore does not affect our estimates of the interaction parameters. 

These values are reasonably consistent with the gas phase Mn value J ~ 0.5eV and with 

n n 

the band calculation [9| which suggests 2A ~ 0.53eV and J ~ 0.65eV (note that in |9j] the 
energy difference for 1 electron with spin up or spin down is 2J with the fit J ~ 1.3eV, 
while from Eq([T]) the difference is 4 J due to \Sc\ = 3/2, therefore one has to divide the J in 
9j by 2 to compare the J fitted here). To summarize, the data and other information are 
consistent with the estimates U = 2.3 ± 0.3eV, 2A ~ J ~ 0.5eV. 
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The analysis of the optics given above was based on the atomic hmit. We have employed 
the semiclassical solver devised by Okamoto et a l |l5 | to solve the model specified by Eq([T]), 



using the tight binding dispersion from Ref 



SJn 



implied by the band theory calculations of ReffQ]. In Fig([T]) we show our calculated results 



with the bandwidth W 



6t 



3eV 



ll| for the change in conductivity across the Neel transition. The results are presented 
in a form which allows direct comparison to the experimental results of Kovaleva et al ^. 
The excellent agreement of energy scales and reasonable agreement of form and magnitude 
confirm the validity of the 
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FIG. 1: (Color online) The change in optical conductivity across the Neel transition. The states at 
temperature T = 0.02t and T = 0.07t are A-type antiferromagnetic and paramagnetic respectively 
111 ]. The plot is designed to be directly compared to the measurement reported in Ref^^], Fig(2). 



We have used the semiclassical approximation to calculate properties for different temper- 
atures, bandwidths and interaction strengths. Representative results are shown in Fig([2]). 
Reducing the GdFeOs distortion is equivalently increasing the bandwidth W. Forcing the 
solution to be orbitally disordered, we found that in the paramagnetic phase the system 
displays a gap at low temperature T = 0.04t for W = Wq (solid) and l.lIFo (dashed) in- 
dicating a Mott insulating state. When the bandwidth is W = I.SWq (dotted) or larger 
the system becomes metallic. We also consider a ferromagnetic state which has an effective 
bandwidth roughly 1.4 (v^) times larger than that of the paramagnetic phase Q,Q. We 
see that for bandwidth W = Wq, the system is metallic for ferromagnetic, orbitally disor- 
dered phase (heavy dash-double dotted curve in Fig(l2])) which is consistent with the results 
from directly varying the bandwidth toW = 1.4IFo (dash-double dotted curve) in the para- 
magnetic phase. Our calculation thus indicates that LaMnOs is a Mott insulator with the 
local interaction strength very close to but slightly stronger than the critical value of Mott 



T 




Frequency t (0.5eV) 

FIG. 2: (Color online) The optical conductivities for bandwidth ranging from W = Wq = 3eV 
(solid) to W = IAWq (dash-double dot). All curves are computed at temperature T = Q.OAt 
which is roughly 30% the calculated orbital ordering temperature, i.e. roughly of order room 
temperature. The heavy dash-double dotted curve is calculated for a hypothetically ferromagnetic, 
orbitally disordered phase, all others are calculated for the paramagnetic but orbitally disordered 
phase. 

transition, and also indicates that ferromagnetic and the large Jh limits discussed above 
provide a poor representation of the physics. Finally we note that if the electron-lattice is 
"turned off", i.e. setting A = in Eq(l), the material becomes metallic (curve not shown). 

At this point we it is useful to discuss more carefully what is meant by the term "Mott 
insulator" in the multi-orbital LaMnOa context. The Cg manifold has a four- fold local de- 
generacy (2x spin and 2x orbital) so an ordering with at least a 4-site unit cell would be 
required to produce a "Slater" insulator. A priori one could discuss a Mott transition in the 
fully symmetric (4x local degeneracy) case, or in the partially ordered case (2x local degener- 
acy) where only one of the spin and orbital symmetries is broken. However our calculations 
indicated that in the manganites the spin states are strongly split so the low energy physics 
is characterized only by a 2-fold orbital degeneracy (along with the global configuration of 
core spins). The criterion for Mott insulator is that the effective inter-orbital interaction 
within the high-spin, orbitally degenerate manifold is large enough to open a gap at zero 
temperature within the single site DMFT approximation. We find that in the paramagnetic, 
orbitally disordered phase the interaction is slightly larger than the critical value. However 
this interaction is composed of two physically distinct contributions, one from interorbital 
Coulomb interaction and the other from the electron-lattice coupling. Removing the contri- 
bution from the electron-lattice coupling changes the interaction strength from slightly larger 
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than the critical value to slight smaller. While in the literature the term "Mott insulator" 
is used for Coulomb-driven phenomena, we think it is appropriate to use it for insulating 
behavior driven by any local interaction, and therefore we identify LaMnOs as being a Mott 
insulator. 

We emphasize that our findings place LaMnOs close to the edge of the insulating phase 
boundary. For this reason details including the presence and evolution with doping of 
GdFeOa rotation (which changes the bandwidth at the 30% level) as well as the possible 
presence of ferromagnetic order (leading to a 40% increase in bandwidth) become very 
important. Finally we note that the semiclassical method used here places the transition 
of the one orbital Hubbard model at U = Ud- In the one orbital Hubbard model the 
difference between Ud and T = metal-insulator point Uc2 ~ l-^Ud is due to the formation 
of a Kondo resonance, which is expected to be suppressed by the core spin coupling in the 
present model. However, further investigation of this point would be desirable. 
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